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The qubit reveals a qubit-cavity system 
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We show that the coupling between a qubit and a single mode cavity can be estimated from 
the process tomography data for the qubit alone. All the parameters of the coupling Hamiltonian 
between the qubit and the cavity mode can be obtained from observations on the qubit. We also 
show that the variance matrix and the photon number statistics of the single mode cavity can be 
reconstructed. Our results provide an alternate means of obtaining the coupling and reconstructing 
the state of the cavity mode in comparison with the techniques presently used. 
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I. INTRODUCTION 

A single qubit interacting with a single quantized radi¬ 
ation mode is the prototypical problem in quantum elec¬ 
trodynamics. The paradigm of cavity quantum electro¬ 
dynamics (Cavity QED) [1] made this experimentally re¬ 
alizable wherein the interaction between a qubit and an 
optical or microwave cavity with a single mode can be 
studied in a controlled manner. Micro-resonators with 
high Q-factors [2] are widely used to prepare single mode 
cavities. [3] In recent years, prescriptions have been pro¬ 
vided to characterize the cavity and also to prepare cav¬ 
ities in desired initial states [4, 5]. 

A relevant and often unknown parameter in a typi¬ 
cal cavity QED experiment is the coupling between the 
qubit and the cavity. Calculating the coupling from first 
principles [6] using parameters such as the atomic dipole 
moment, field amplitude and the volume of the cavity [7] 
is typically based on several assumptions and idealiza¬ 
tions, necessitating a direct estimation of the coupling in 
a real experiment. This is done [8-11] by observing the 
Rabi oscillations between the normal modes of the sys¬ 
tem [12, 13]. The splitting between the levels is propor¬ 
tional to the coupling at resonance, and it is also related 
to the excitation quanta of the cavity. The experimen¬ 
tally estimated coupling strength may deviate [9] from 
the number predicted by the cavity geometry [7]. 

In this Paper we investigate a method of estimating 
the qubit-cavity coupling from the observed dynamics of 
the qubit alone in the strong coupling regime. A simi¬ 
lar question was previously addressed for the case of two 
coupled qubits [14] as well as for a qubit coupled to an N 
level system [15] which showed that the coupling could be 
obtained from the process tomography data of the qubit 
alone. We also obtain the variance matrix [16] of the ini¬ 
tial quantum state of the cavity. Our recipe can be ex¬ 
tended to obtain the higher order moments of the initial 
state of the cavity as well, leading to a complete quantum 
state tomography of the cavity mode [17-19]. While state 
tomography of the cavity mode has been done [18, 20, 21], 
we provide an alternative approach for doing the same. 
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Compared to the available procedures for estimating the 
coupling [8-11], our prescription does not demand any 
prior knowledge on the excitation quanta of the cavity. 


II. A QUBIT COUPLED TO SINGLE CAVITY 
MODE 

The state of the qubit is described by the Pauli matri¬ 
ces, CTi, CJ 2 , CTa which satisfy the SU(2) Lie algebra 

[fjj, (7j\ = 2iCijj„(7j„ 

The single mode cavity has a characteristic frequency to 
and it is described by the bosonic operators b and 
satisfying the usual commutation relation, 

[5,6t] = i 

First, we consider the strongly coupled case where the 
atom and the cavity mode are near resonance. In this 
limit, the strong interaction between the two is gov¬ 
erned by the Jaynes-Cummings [22] Hamiltonian. In 
Section. IV, we also consider the case when the qubit 
and the cavity are weakly interacting [23, 24], which is 
the off-resonant or dispersive limit of cavity QED. The 
Jaynes-Cummings Hamiltonian is given by 

H = ‘^<73+uj(^b'^b+^ + g{ba+ + b'^a_), (1) 

where a± = (cti ± %a-2)l‘2.. Our aim is to estimate the 
parameters a, a; and g appearing in H. We have used 
atomic units (h = 1) in the above expression and the 
rest of the Paper. 

To obtain the parameters, the qubit is initialized in the 
following three states, 

+ fc = l,2,3 

At the moment, we do not worry about the specific ini¬ 
tial state for the cavity, even though prescriptions are 
available to initialize the cavity to a handful of possible 
states [4, 5]. 
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The state of the qubit, in the Schrodinger picture, at 
any later time t can be written as 


Similarly the other two first derivatives can be evaluated 
as 


pf ^ + ci^\t)a2 + 4^\t)a3] 

(k) 

Here, pl is the density operator of the qubit obtained 
by tracing out the cavity, for different initial states of the 
qubit 

The dynamics of the qubit is specified by the nine 
quantities which we assume are obtained from a 

process tomography experiment corresponding to a du¬ 
ration t on the qubit. Notice that are obtained as 

the expectation values of Pauli operators with respect to 
the specific initial state as 

cfV) = (0-*)^^ = = tr[77f^cri (g) Ic] (2) 


Here Ic is the identity operator acting on the cavity. 
Since full or partial traces of the identity operator with 
respect to any state of the cavity is equal to unity, we 
avoid writing Ic explicitly in the following discussion. 
The second trace in Eq. (2) is evaluated over the den- 
sity operator for the combined system rjl . Switching 
to the Heisenberg picture, where the time dependence is 
associated with the operators, we notice that 

=tr[?7^''Vi(t)] 

(k) 

In the above expression, tjq is the initial state of the 
qubit and the cavity taken together. Since we have pre¬ 
pared the qubit in a specific initial state, we know that 
the combined initial state is separable. 


(k) (k) „ (c) 

Vo = Po ‘^Po 


( 3 ) 


Here p^^ is any arbitrary quantum state of the single 
mode cavity which can be repeatedly prepared. Later on 
we show that our treatment can be used to obtain partial 
information about plf\ 

In the Heisenberg picture, the time derivative of 
is. 


d” 





k 


= tr 




( 4 ) 


(T 2 (t) = acri(t) - g{h + h^)a 3 {t) 

= gib + d'f)cr 2 (t) - igib - b'^)(Tiit) 


Using these relations we calculate the first derivative of 
(t) defined in Eqn. (4) at t = 0. The initial state of the 
qubit-cavity system is separable as in Eqn. (3) and the 
trace appearing in Eqn. (4) can be decoupled and com¬ 
puted as product of traces over the qubit and the cavity 
mode separately. We obtain the following relations, 

Ci^^(O) = —a 
cP^(O) = ig{b-b^) 

4^^(0) =+ (5) 

We notice that the tensor is anti-symmetric in i 

and k. From Eq. (5) we obtain, 


a = -cf4o)> 


(6 + &t) = -- 


^9 

a(3) 


( 0 ) 


( 6 ) 


The first derivatives have fixed one Hamiltonian parame¬ 
ter, namely the free evolution parameter, a, of the qubit. 
This is an expected result since we are looking at the 
dynamics of the qubit system alone. 

In the same spirit, we also ask for the second deriva¬ 
tives of The second derivative of cri(t) is 


(j'lit) = i[H,di{t)\ 

= -a^cri(t) -h agib + b^)cF3{t) - g‘^ii{b - b^)f(Ti{t) 
+gu!ib + h^)a3it)-‘rig^ib - 6^ )cr2(t) 


The second derivatives of other Pauli matrices are also 
computed in a straight forward manner. 


d 2 (i) = -a^U 2 it) + iagib - b^)(J 3 {t) - g'^{h + b^f'(J 2 {t) 
+igLoih-b^)u3it) + ig'^ib -h^ )<^iit), 


The derivatives can be computed from the time series 
data available from process tomography experiments on 
the qubit which makes the scheme experimentally realiz¬ 
able. In what follows, we will show that the time deriva¬ 
tives of (t) contain extractable information regarding 
the Hamiltonian parameters and the initial state of the 
cavity. 

We start with the first derivative of cti (t), 

diit) = = f[H, e*^Vi(0)e“*^‘] 

= fe*^*[iL,ai(0)]e-*"‘ 

= -a<T 2 it) + igib - b'^)a 3 it). 


and 

d’3it) = gaib + b'^)aiit) - g'^ib + b'^fa3it) 

— igujib — b^)a2it) + igaib — b^)cr2it) 

- g^iiib - h'^)fa3it) - gujih d- 6^)cri(t) - ^g"^ 

We employ Eq. (4) again to compute the second deriva¬ 
tives of From tT'i(t) evaluated at t = 0, we obtain 

the following relations. 

ci«(0) = -a2-5"((*(6-6t))2), 

Ci(2)(0) = ig'^{b - 
ci^^^(O) = agib + 6'^) -|- gujib + 6'^). 


( 7 ) 
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Another set of equations are obtained from (T 2 {t) as, 
c,i^)(0) = -a^-g^(b + b^f), 

= iag{b — b^) + iguj{b — 6'^). (8) 

The remaining set of relations are obtained from d'sit): 

C 3 ^^^( 0 ) = ga{b + b'^) - guj{b + b'^) - 2g'^, 

= —iguj{b — 6'^) + iga{b — 5'^) — 2g^, 
cs^^Ho) = + b^r) - g^{i{h - b^)f) - 2g\ (9) 

Though most of the relations are redundant here, they 
are useful in checking the consistency of the prescription 
we provide. Partial information regarding any of the pa¬ 
rameters will also help us to decide the resolution of the 
time series data needed for a demanded accuracy on pa¬ 
rameter estimation. 

The second derivatives and evaluated 

at t = 0 can be used to compute the coupling g as 

9=^1 [c'i(')(0) + C2(')(0) - C3(')(0) + 2a2]. (10) 

We obtain the expectation values (6±6l) by substituting 
Eq. (10) into Eq. (6). These expectation values are pro¬ 
portional to the mean values of the x and p quadratures 
of the cavity mode with 

x=-^ib + b^) and p=^(6-5^). 

If the cavity mode is in a Gaussian state, then the state 
is fully characterised by the variance matrix defined as 

( (A2) {Ux,P}) 

V= 1 _ 2 ^ 

(P2) 

with X = X — (x) and P = p — (p). By inverting a subset 
of equations from Eqs. (7), (8) and (9) above and using 
Eq. (6) we obtain the elements of the variance matrix of 
the cavity mode as, 

V.. = (X^) = - C2(")(0) - (C2(^)(0))^] 

Vpp = (P^) = 2^ [ - - {cP^\0)r] 

III. NUMERICAL EXAMPLE 

We have simulated a simple numerical example for the 
qubit-cavity system in the absence of real experimental 
data. For this purpose, we assumed the Hamiltonian pa¬ 
rameters a = 1, w = 1 and g = 1. We also assumed 


that the cavity is initially prepared in a coherent state 
with (TV) = 1. The time series data for different initial 
preparation of the qubit state and were 

constructed from the exact Hamiltonian evolution using 
MATLAB. For numerical evaluations, the number of lev¬ 
els in the quantum oscillator were reduced to four hun¬ 
dred. The data points were taken at discrete points to 
mimic a real experiment, by evaluating the reduced den¬ 
sity matrix for the qubit at those points. An additional 
noise could also be added to the data points, but no such 
noise has been added to the time series data we used 
here. 

With the artificial time series data we constructed, we 
computed the first derivatives Ci^^^(O) and the second 
derivatives Ci^^'^(O) at t = 0. The derivatives evaluated 
at i = 0 are only approximates to the true value of the 
derivatives. They were computed using a finite but small 
step size 6, which cannot be arbitrary small in a real 
experiment. Typically, the smallest step size is the time 
taken by the qubit to pass through the cavity after its 
initial preparation. 

We then used the relations obtained in the previous 
sections to reconstruct the Hamiltonian parameters as 
well as the variance matrix of the cavity mode. The re¬ 
constructed Hamiltonian parameters are shown in Fig. 1 
and the variance matrix elements are shown in Fig. 2. 
We observe that the maximum relative error is less than 
two percentage, for a step size of 0.01. Notice that the 
natural timescale in the system, set by the inverse of the 
Hamiltonian parameters is also equal to one, since we 
have set all the Hamiltonian parameters equal to one. 



FIG. 1. Reconstructed Hamiltonian parameters as a function 
of the step size S used to compute the derivate. The true 
value assumed is unity. 
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Eq. (12) can be inverted to obtain 





Var(7V) = {N^) - {Nf 
4g4 




4|^+4af)W 


21 


9 , ry 9 

a5+"“a 


-(N)- 


This procedure of evaluating the higher derivatives of the 
time series data near t = 0 can be extended to obtain 
the desired higher order moments of the photon number 
distribution of the initial state of the cavity. 


FIG. 2. Reconstructed variance matrix elements for the co¬ 
herent state with photon number one, for different step sizes 

S. 

IV. DISPERSIVE LIMIT 

Another regime of interest in Cavity QED, is the dis¬ 
persive (off-resonant) limit. The system is then described 
by the following Hamiltonian [23]: 

2 

H = ujiN + 1/2) + “as + + l/2)a3 

= a;(V + l/2) + ^a3 (11) 

Here A is the atom cavity detuning given by A = a —w 
and [H, A] = 0. Note that the a and g appearing in this 
Hamiltonian is the same as in Eq.(l). This weak coupling 
between the qubit and the cavity has been used to resolve 
the number states of the cavity [24]. Here we show that 
our method can be used to probe the photon number 
statistics of the cavity by revealing the mean and the 
variance of this distribution. Here we assume that the 
parameters a and g have already been estimated using 
the methods outlined above. 

We notice that CTi(t) = —Aa 2 {t) and & 2 {t) = Aai{t). 
Similarly, the second derivatives are evaluated to obtain 
di(t) = —A‘^ai{t) and 0 - 2(0 = —A^a- 2 (t). These relations 
give us the following expressions: 

cW(0) = (i) = a + 2^(^(V) + i) 

-cW(0) = (iO = 4^(iV2) + ( 4 ^ + 4a^) (N) 


V. CONCLUSION 

In this Paper, we provide a scheme for estimating the 
Hamiltonian parameters and for obtaining the variance 
matrix for the quantum state of the cavity by observing 
the dynamics of a qubit in a cavity-QED system. In ad¬ 
dition to obtaining the variance matrix [16] which com¬ 
pletely characterizes Gaussian states of the cavity, our 
recipe can also be extended to obtain the higher order 
moments of the initial state of the cavity. The variance 
matrix of the initial state of the cavity mode is useful for 
designing quantum information processing protocols us¬ 
ing continuous variable quantum systems. We also stud¬ 
ied the dispersive limit, and we have provided an alterna¬ 
tive prescription to resolve the photon number states of 
the cavity [24] by measuring the qubit. The calculations 
presented here can be repeated for a situation when there 
are N cavity modes interacting with a qubit, with char¬ 
acteristic frequencies uji and couplings gi . Our treatment 
of evaluating the time derivatives of the process tomog¬ 
raphy data up to second order in this case would give 
us the sum of squares of the couplings 
further details about the distribution of the couplings of 
the N modes to a qubit. 
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